close all; clear all; clc; 


psi1=0.3471;psi2=2.44;

psi1b=0.3465;psi2b=2.08;

a = load('agrid2.txt');
a = exp(a);

udef = max(0,psi1+psi2*log(a));
udefb= max(0,psi1b+psi2b*log(a));

figure;
plot(a,udef,'r',a,udefb,'b');
